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An effective Hamiltonian, describing quantum tunneling in ferrimagnetic nanoparticles which 
includes interactions between the electronic spins of nanoparticle and microscopic environmental 
spins (like nuclear spins or paramagnetic impurities), is obtained. Two limiting cases, describing 
tunneling in antiferromagnetic and ferromagnetic regimes are considered, and criterion for the tran- 
ce ■ sition between the two regimes is found. The validity of analytic results is verified by the exact 
diagonalization method. 

G^ ■ 



I. INTRODUCTION 



Q 

The quantum tunneling of magnetization (or Neel vector) has been attracting a lot of attention during the last years 
from both experimental and theoretical point of view since this quantum effect might be seen on macroscopic scale for 
the total number of spins of order 10** and more. Electronic spins in the magnetic grain are coupled to each other by the 
' exchange interaction, forming either a "giant spin" in ferromagnetic case or a "giant Neel vector" in antiferromagnetic 
; case, and can tunnel coherently between the two degenerate states separated by the magnetic anizotropy barrier. The 
theory of coherent tunneling of the magnetization vector has been developed in many papers [1] (for a review see [2] 
] and refs therein). As to the antiferromagnetic particle, the tunneling effect in this case may be even stronger than 
C\| ■ in the ferromagnetic one [4] (for a review of modern state of the theory see [3-9]), and of the experiment - [10].) It 
' turns out that in both cases interactions with a spin bath (like nuclear spins or paramagnetic impurities) play a very 
important role [11,13]. These interactions may be so strong that they can completely suppress the quantum coherence 
[12]. For the theoretical analysis of the tunneling dynamics one usually introduces an effective Hamiltonian, H^ff, 
operating only in the low-energy subspace of the combined system - "giant spin" (or "giant Neel vector" ) plus the 
spin bath. This procedure involves truncation of the original problem to the subspace of the two lowest energy levels 
5^ , corresponding to opposite spatial orientations of the order parameter, i.e., all transitions from the lowest doublet to 
states with exitation energy E > fto, where fto characterizes the small oscillation frequency of the order parameter 
^ near the classical minima, are ignored in H^ff- Earlier -ffe// for the "giant spin" was derived by the instanton 
O technique in [15]. In this paper I obtain an analogous Hamiltonian for the case oi ferrimagnetic/ antiferromagnetic 
^ nanoparticles. To explore the limits of validity of the analytical expressions I use the exact diagonalization technique. 
> 

^ ■ II. HAMILTONIAN AND MODEL 

; A. STARTING MODEL 

Taking advantage of the known fact that Heisenberg model with antiferromagnetic ordering can be reduced in the 
long-wave limit to the non-linear 0(3) cr-model (both in ID and 2D cases) [3,5,7,9,16], we introduce an effective 
two-sublattice model with strong exchange interaction between sublattices, the Lagrangian being analogous to that 
of 0(3) cr-model (see, for example, [3,4,8,18]). In addition to the just mentioned works we include coupling between 
the sublattice spins Si and 52 and environmental spins (| ai \= i) in a weak magnetic field Ho- Since coupling to a 
single environmental spin is typically rather small and can be treated perturbatively, it is sufficient, without loss of 
generality, to consider here just one such spin coupled, say, to the sublattice 5*1. In this paper we concentrate on the 
contact hyperfine interaction, assuming that it is a nuclear spin. Then, choosing the easy axis/easy plane magnetic 
anizotropy, we write the starting Hamiltonian as follows: 

H = JSi-S2 + U{Si) + U{S2) + ■ a, U{S) = -K«Sl + K^Sl ~ jeHoS, (1) 

where lOo is the hyperfine frequency, K^^ > and K± > are magnetic anizotropy constants, and we assume that 

J, K± ^ 7e'5'i,2 I Ho |<C 2i4r||5'i 2, ^ 2ii'||5f 2- (We do not confine ourselves hy J :s> K±, although in practice 
this limit realizes very often). Assuming Si + S2 ^ 1, wc will use the semiclassical description. 
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In the low-energy dynamics the dominant contribution to the transition amphtude between the two semiclassical 
energy minima comes from trajectories with almost antiparallel Si and 82- That is Si and S2 tunnel simultaneously, 
with only small fluctuations about the comon axis, and the tunneling is described by the Neel vector N = Si — S2 
with almost constant modulus. As we assume in the general case that Si ^5*2, an excess spin S = Si — S2 will tunnel 
together with the Neel vector (ferrimagnetic case). 

Let us introduce the spherical angles 9i, 62, <t>i, <f>2 to describe the orientation of 5*1 and 5*2. We then write 
62 = n — 61 — eg, <j)2 = 7r + 0i +£0 with jegl, le^l <S 1 (see [4]). Also, we choose as the basis set the states characterized 
by the opposite orientations of Neel vector (along "easy axis" z), | JJ.), | ff). Then we can write the transition amplitude 
as the path integral: 

r„/3(i)=^ I?{(?l,02,</'l,02}exp|-^ dT[Co{T)+5CH{T)+5C„{T)]^ , (2) 

where Co{t) - Lagrangian of the non- interacting problem, 5Ch{t) describing the weak magnetic field contribution 
and 5C„{t) describing the interactions with the nuclear spins (all the Lagrangians are in the Euclidean form and 
a,/3=|4),l^)). 

In our problem there are two time scales. Q,^ ^ corresponds to the transition time between the two semiclassical 
minima. The second time scale is defined by the value of the tunneling amplitude Aq (the tunneling splitting is 
equal to 4Ao cos -kS) . The relation between the transition amplitude (|^) and the non-diagonal part of the effective 
Hamiltonian is given by [15] (for the real-time axis): 

(i7f^y(f))^^ = ^r°^(t) ; (a^/3, l]-i«t«A-i). (3) 



B. SEMICLASSICAL SOLUTION 

Let us find first an extremal trajectory for Oi{t) and 0i(t) for the non-interacting problem when \Ho\ = and 
LJo = 0. The Lagrangian Co has the following form: 

Co = JS"! 52 (sin 6*1 sin 6*2 cos(0i - (/12) + cos6ii cos6'2 + I) + C% + C% , 

C°g^ ~iS(t)e sine + K\iS^ sin^ 9 + K^S^ sin^ e sin^ <j). (4) 

Since y-axis is a "hard" one, (j) will weakly fluctuate around (j) = or (j) = tt. These values of (j) correspond to 
two possible trajectories connecting two semiclassical minima: clockwise and anticlockwise. Performing the Gaussian 
integration over three "fast" variables (0, ee, e^) we get (omitting all the terms which do not contribute to the equation 
of motion): 

Co{e) = + ^.^' Sin' 9 + K«sm^e, ~z = Ki 



2 2J " ' J + 2K^Sl' 

I ^SSIK^ ^ Si ^ 2SIKI 
2S z{J + 2K^Sl) J + 2K^Sl i{J + 2K^SlY' 

where = K^ \\{Si + 5|) and J = JS1S2 (below we will ignore indexes at 9 and (j)). This Lagrangian results in 

the equation of motion 9 ~ ^lo sin9, the solution of which is: 

sin 6I(t) = l/cosh(rJoT), Vto ^ {2k\\/ M^''^ , (6) 

Substituting the extremal trajectory into Co and integrating over r, we get for the Euclidean action: 

A'l ^Ao + iwS, Ao = 4A'|| /ilo, (7) 

where 77 = ± corresponds to the motion between two semiclassical minima in clockwise or anticlockwise direction and 
rjirS is the Kramers-Haldane phase. 



C. MAGNETIC FIELD 
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Now we introduce a weak magnetic field {"fe = 1)^ 



5Ch = SCf^ + 5Cl^ , SC^i = -S{H^ sin 6* cos </) + Ry sin 9 sin 4> + cos 61). (8) 

If we integrate Co + 6Ch over ee,e^, with the excess spin 5 equal to zero and take the limit J ^ K±, then we will 
get the well-known expression for the Lagrangian of Andreev-Marchenko [17] without the gradient term and the term 
describing weak ferromagnetism of Dzyaloshinskii. But we will proceed in the general form. After integrating (Jq) over 
fast variables eg, and </> (keeping only linear in magnetic field terms) we substitute the extremal trajectory (m and 
further integrate over r to get a correction to the action (Q): 



" J + 2K^Sl\ i{J + 2K^Sl)] no \ J z(J + 2A'i5|)2 

r^^S \ Sno \ ynSnoSl \^ , . 4j^±^2 ^ 1 



\ 2~z 2~z{J + 2K^Sl) I no 



D. NUCLEAR SPINS 

Let us include now coupling between the electronic spins of the nanoparticle and nuclear spins. The correction to 
the Lagrangian Co from hypcrfinc interaction has a form: 

5Ca = (^x sin 9i cos + dy sin sin (pi + dz cos 6*1^ , (10) 

where Ui {i = x,y,z) are Pauli matrixes. Performing the same operations as in the case of magnetic field we get for 
the correction to the action: 

xAn m^of. , . SlnoKx_ . \ .muJoS , 

oA ' = ax + I — ^ cr„ + « cr„ . 11) 



E. EFFECTIVE HAMILTONIAN 

Since n~^ <^ ^o^ ^ we can write the non-diagonal part of the effective Hamiltoniah in the form: 

= \ {r^t^^it) + H.c] , r^^(t) = ^tAo ^ exp {-A"} , 
where f_ is the lowering operator in the space of Pauli matrixes and the transition amplitude is (see also [14]): 



(12) 



Ao = ^J-Aoexp{-Ao} 



(13) 



We start constructing the effective Hamiltonian from the antiferromagnetic case, when 5*1 — S2 and only Neel 
vector undergoes tunneling. Keeping only linear in Ho and lOo terms in the action (for the simplicity we set = 0) 
we get: 



^"(0) = Ao{0) - 17^^(0) + r/a(0) ( <t, + iX{0)ay J , ri„(0) = 282^ K\\{J + Ki_) 



Ao{Q) = 452A(0), V(0) 



«(0) 



J I 



J + K\' 2rJo(0)' 
Substituting (|l|) into (|l|), we obtain the non-diagonal part of the effective Hamiltonian: 



ijf.?(0) = 2Ao(0)f_ cos 



V'(O) -l-a(O) i&x - A(0)(T. 



H.c. 



(14) 



(15) 
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We do not present here the diagonal (static) part. The general way of calculating this value one can find in [15]. 

It is easy to verify that in opposite limiting case, Si ^5*2, we will get the effective Hamiltonian describing tunneling 
of the ferromagnetic-grain magnetization: 



H^jf (Si) ^ 2 AoiSi)T^ cos 



n^iSi) = 2SiJK»K^, a{Si) 



nSi - ip{Si) - a{Si) ( iax - KSi)^y 

ttH. 



H.c. , 



2no{Si) 



All 



2K 



" AoiSi)^2Si\{Si), KSi) = \Iy^, (16) 



where is understood as the grain's " giant spin" . 

Finally, in the case of arbitrary value of S (ferrimagnetic) the Neel vector tunnels together with the excess spin. In 
the limit of J 2K±S2 the effective Hamiltonian has the following form (as before = 0): 



HeffiS) = 2Ao(S')f_ cos 
Ao{S) = 



nS — tp{S) — a{S) yi<Tx ~ A(S')(7y 



H.c. 



52 + 51^ + 553 



J 



TTAo{S)Hy _ TTLdp 

2no{s) ' " 2no{S) 



J J 



, A(5) = ^f^ + ^Y 



AKuK±J 



(17) 



From the analysis of these three Hamiltonians the criterion (on the value of the excess spin), corresponding to the 
switching between the ferromagnetic and the antiferroniagnetic regimes, follows directly: 



5»52 



5«52i 



2Ki_ 



J + 2K^Sl 

I 2WI 
J + 2K±S^ 



ferromagnetic regime 



— anti ferromagnetic regime. 



(18) 



III. EXACT DIAGONALIZATION 

To verify the obtained expressions we use the exact diagonalization method in application to the starting Hamilto- 
nian (|l|) with Hx — 0. Because of the limited space, we present the numerical results only for the antiferromagnctic 
case (5i = 52). The exact diagonalization procedure is described in detail in [15], and therefore we define here only 



the values that are subject to the analysis. Expanding (15) into the series up to the linear in lUq terms we get 



H^ff w 2Ao(0)|^cosi/>(0) • T^ao - a(0)sini/' • Tya^ + a(0)A(0) sin?/; • f^ayj. (19) 
An effective Hamiltonian defined by the exact diagonalization has the form: 

+ Hd , (20) 

where Hd is diagonal in the f^-representation. The coefficients Cg are presented in Figures l.a, l.b and l.c (dashed 
lines) together with the corresponding analytical values from (nS) (solid lines) as functions of magnetic field at 
J = 40, K± = 20,A'[| = l,Wo = 0.2, 5i = 52 = 10,cr = 1/2. In Fig. l.d the coefficient a extracted from Cy^ 
(dashed line) is shown together with its analytical analog from ( p^ ) (solid line) as a function of sublattice spin at 
J = 60, K± — 20, = l,CiJo = 0.2, a = 1/2. This coefficient describes the interaction with the spin bath and, 
hence is the most interesting for the present investigation. As one can see from all the Figures, the coincidence of 
the analytical results with the numerical ones is quite satisfactory. Some differences between the two values for a at 
small 52 are explained by the fact that semiclassical description is far from being valid in this region. 
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IV. CONCLUSIONS 



Though the detailed analysis of all the factors effecting the tunneling process is beyond the scope of this paper, it is 
possible to make some qualitative conclusions (even without rigorous calculations with obtained Hamiltonians) on the 
spin dynamics of the investigated systems: (i) the quantum tunneling effect is really more strong in antiferromagnetic 
nanoparticles than in ferromagnetic ones that can be seen even from the fact that Ao{0) < Ao{Si); (ii) the presence 
of the excess spin strongly effects the quantum tunneling process due to the mere fact that it renormalizes Ao and 
which both define the tunneling splitting and the bounce frequency; (iii) interaction with the environmental spins is 
extremely important since it can drastically change the tunneling picture down to the full suppression of it, as in the 
case of tunneling of half-integer excess spin or half- integer " giant spin" . 
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FIGURE CAPTIONS 



Figure 1 (a)-(c) The numerical cocfBcients Cij (dashed hne) in comparison with the analytical expressions from 
(H) (solid line) as functions of magnetic field Hyat J = 40, K± = 20, = 1,lu^ ^ 0.2, Si = S2 = 10, ct = 1/2. (d) 
The numerical (dashed line) and the analytical ( |l9| ) (solid line) dependences of the coefficient a on the value of the 
sublattice spin at J = 40, Kl. = 20,K\\ = l,ujo = 0.2, a = 1/2. 
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